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Abstract 

The extrapolation of the decay amplitudes of the pseudoscalar mesons 
into two photons from the soft meson limit where it is obtained from 
the axial-anomaly to the mass shell involves the contribution of the 0~ 
continuum. These chiral symmetry breaking corrections turn out to be 
large. The effects of these corrections on the calculated n° decay rate, on 
the values of the singlet-octet mixing angle and on the ratios j 2 - and j 2 - are 
discussed. The implications for the transition form factors 77* — > ir°, 
•q, 7/ are also evaluated and confronted with the available experimental 
data. 



1 Introduction 

Currrent Algebra and PCAC (partial conservation of the axial-vector current) 
together with the ABJ axial-anomaly [1] expresses the two- photon decay rate 
of the pseudoscalar mesons (P = 7r° , 77, 77') in terms of the coupling constants 
f^o , fo and fg and of the singlet-octet mixing angle 9 

(0\A a \ 7T°(p)) = -2ifop a 

(0 I j „(p)) = 2if 8 . co S 9. Pa , (0 I I V (p)) = -2</ . sinf?.p Q 

(0\AW \T,'(p))=2if 8 .sm0.p a , (0\AW\r!'(p))=2if .cos8.p a (1.1) 
with 

~ /,+ = 92.4MeV 

and where the axial- vector currents are given in terms of the quark fields. 
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Unlike /jr, /o and /gare not directly related to any physical process. 

SU(3) breaking enters through singlet-octet mixing and through the devi- 
ation of the ratio y 2 - from unity.Another symmetry breaking effect originates 
from SU(3)xSU(3) breaking: the decay rate is obtained in the soft meson limit 
from the ABJ anomaly [1] and the extrapolation to the mass-shell involves cor- 
rections 0(m 2 P ) which are expected to be small for the ir° but which are not 
necessarily so for the rj and r\ These corrections to the PC AC limit, which arise 
from the CP continuum are similar to the corrections to the Goldberger-Treiman 
relation [2]. Attempts to estimate these corrections have been undertaken [3], 
[4], [5]. In [3] and [4] only the contribution of the high energy part of the 0~ 
spectrum was considered and in [5] only the low energy part was used. These 
three calculations are moreover heavily model dependant. 

It is the purpose of the present work to provide an estimate of the PCAC 
corrections which enter in the evaluation of the two photon decay rates of the 
pseudoscalars under the sole assumptions that the V-A-A vertex is given by the 
quark triangle graph in the deep euclidean region and that the main contribution 
to SU(3)xSU(3) breaking arises from the energy interval 

lGcV 2 < s < 2GeV 2 

of the 0~ continuum. 

In section 2 we give details of the calculation and the resulting constraints on 
the values of /§, /o and the mixing angle 9 are discussed and finally an evaluation 
of the implications of our results on the transition form factors 7*7 — > P are 
presented in section 3 and confronted with the available experimental data. 

2 The chiral symmetry breaking corrections 

Consider the three-point function 
T a „ v {p,qi,q 2 ) = 2m. J J dxdy exp(-iq lX + ipy)(0 \ TA a (y)V„{x)V v (0) | 0) 



e^x a qUZPa-F{p 2 ,qlql) + ... (2.1) 



7r (p 2 — m 2 ) 

where the dots represent other tensor structures and where V^„ denotes the 
electromagnetic current 



2 1- 1 

^ = 3^7^ - gtffyd - gS7 M s (2.2) 
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The pion pole contribution having been isolated in expression (2.1). 

In the soft it limit and with both photons on mass-shell, the axial-anomaly 
[1] yields 

F(0,0,0) = 1 (2.3) 

The rate of the decay ir° — > 2j provides a measurement of 

F(n4,0,0) = l + A^A* = 0(m*) (2.4) 

Equation (2.1) is used to define the off-shell symmetric amplitude 

F(s,t)=F(s=p 2 ,t = q 2 = q 2 2 ) (2.5) 

In the deep euclidean region F(s,t) = F QCD (s,i), F QCD (s,i) includes per- 
turbative and non-perturbative contributions 

F^ CD (a,t) = F p (a,t)+F np (a,t) (2.6) 

The perturbative part F p (s, t) is obtained from the quark triangle graph, the 
contribution of which takes a particularly simple form in the symmetric case [6] 

F p (s, t) = -2.(, - ml). J f Q d X dy. {y{ j;y 2 t)+ty * = 1 - « (2.7) 

The non-perturbative part F np (s,t) arises from the contribution of the vacuum 
condensates 



where 

7r 3 ,a s G 2 . , 64 . . , .,. , , . . 

, b 2 ^~-.^.{a s (qq) 2 ) , 6 3 = -&i (2.9) 

For t fixed, F(s, t) is an analytic function of the complex variable s with a cut 
on the positive real axis which extends from s = 9m 2 to oo [7] . 

A dispersion relation between F(s = m 2 ,t) and F(0,t) is obtained from the 
integral 

hL^<r FM ,2 ' 10) 

where c is the closed contour in the complex plane consisting of a circle of large 
radius R and two straight lines lying immediatly above and immediatly below 
the cut .Cauchy's theorem yields then 

F(^ ( ) = F(0.0 + g^ I ^.F<M) (2.11) 



3 



The integral in the equation above consists of two parts : an integral of 
the discontinuity of F(s,t) over the cut , which provides the main contribution 
and an integral over the circle of radius R where F ( ^ CD (s,t) provides a good 
approximation to F(s, t) except possibly in the vicinity of the positive real axis. 
Little is known about the integrand over the cut which consists of the contri- 
bution of the axial- vector (1+) and pseudoscalar (0~) intermediate states. The 
contribution of the low energy region(the three pion states) is suppressed by 
loop factors and amounts to little [8]. We expect the major part of the contri- 
bution to originate from the al(1260) and 7r'(1300) bumps i.e from the range 
lGeV 2 < s < 2GeV 2 .In order to eliminate this contribution, we shall consider, 
instead of integral (2.10), the following modified integral 



i-ira-<;-°»-°" s >- F <»' ( > < 2 ' 12 > 



In the equation above . the coefficients a and a\ will be chosen so as 
to annihilate the integrand at m\ = rr? al = 1.56GeV 2 and at m 2 , = m 2 , = 
1.70GeV 2 , i.e 

ao = — 5- H j , ai = 5 — 5- (2-13) 

This choice reduces the integrand to only a few percent of it's initial value 
over the interval lGeV 2 < s < 2GcV 2 . It will also considerably reduce the 
contribution to the integral over the circle near the positive real axis. The con- 
tribution of the continuum having thus been drastically reduced, we neglect it. 
Cauchy's theorem now yields 

F(m 2 , t) w F(0, t) + a .m 2 .F(m 2 , t) + a x .m\ .F(m\ , t) 

1 1 .(--a Q - ai . S ).F(s,t) (2.14) 



27ri J (s — m 2 ) s 

The integral being now carried over the circle of radius R. 

In order to obtain F(s, 0) we proceed in a similar fashion. F(s, t) is an analytic 
function of the complex variable t except for a cut on the positive real axis 
extending from t = 4m 2 to 00. In the low t region F(s,t) is dominated by the 
p — u double and single poles, 

VI f\ C l( S ) , C 2(g) ro ,,x 



Consider now the integral 

{t-m 2 p f 



2m J c , 



t 



.F(s,t) (2.16) 



where c' is a contour similar to c in the complex t plane, the double and single 
vector meson poles having been removed and for an appropriate choice of R' 
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the major contribution to the integral (2.16) comes from the integral over the 
circle so that 



1 1 f it — m 2 ) 2 

Fi » A]) ~^-™L*—t jL - Fi "- i) 



It follows from equations (2.14) and (2.17) that 

F(m£,0).(l - a .ml - Oi.m^) 

(2.18) 

The integral above being carried over the circles of large radii R and R 1 we have 
replaced F(s,t) by F QCD (s,i) in the integrand. Inserting eqs (2.6) - (2.8) in eq. 
(2.18) gives 

F{m 2 , 0).(1 - a .m^ - ai.m*) 
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TYl 

+2.-|.oi.63 (2-19) 

The non-perturbative contributions to equation (2.19) turn out to be negligible. 

The decays rj — ► 27 and 77' — > 27 are treated in a similar fashion. The 
axial-vector currents which project on the r\ and rf states are, respectively, 

4(8) AO) AS) AO) 

Al = R^. cos0 sine) , Ai = (^-.sme+^.cos0) (2.20) 

.18 JO .18 JO 

These currents we use in the definition of the 3-point function equation (2.1). 
The results for the two photon decay rates of the pseudoscalars are then: 

^°-^'Mf/ i+ ^ 2l t <2 - 2i) 
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W- 2 ,> = ^.( 1 + A,,^ + 2v^ (, 2 3, 

■f 

With the notation F ,8 = j 2 ^ and 

F(m^0,0) = (l + A 7r ) 

F(m^0,0) = -^(l + A JJ ) , F(m2„0,0) = ^(l + A,0 (2.24) 

A T is obtained from equation (2.19) The same equation with the negligible non- 
perturbative part omitted yields A,, and A^/when m^ is replaced by m^/and 
when the strange quark mass is neglected. 

R', the duality radius in the p- meson channel is usually taken to be R' « 
1.5GeV 2 in the litterature [9] The value of R should be large enough to include 
the contribution of the pseudoscalar excitations but not too large in order not to 
invalidate the approximations at hand. Moreover Ap should be stable against 
small variations in i?..For the parameters m\ and m 2 we take 

ro? = m|| l(1260) = 1.56G6V 2 , ml = m 2 (1300) = 1.70GeV 2 (2.25) 

as discussed previously . 

In the 7] and rj channels two pseudoscalar excitations, 77(1295) and ?7(1440) 
as well as two axial- vectors, /i(1285) and /i(1420) dominate the 0~ and 1+ 
continua and couple to both r\ and 77' mesons with unknown strengths. Because 
77(1295) and /i(1285) on the one hand and 77(1440) and /i(1420) on the other 
hand are practically degenerate in mass it suffices to take 

m\ = 1.66GcV 2 , m 2 2 = 2.04 GeV 2 (2.26) 

for both 77 and r}' . 

With these choices expression (2.19) passes through a maximum for R m 
2GeV 2 in the tt channel and R ^ 2.25GeV 2 in the 77,77' channels., values which 
we adopt because of the stability criteria stated above. Numerically then, we 
find 

A,, - .047 , A,, - .77 , A,- - 6.0 (2.27) 

The relative theoretical error on the numbers above is estimated to be of the 
order of a s (2.5 GcV 2 )/tt=12% 

The chiral symmetry breaking corrections are thus seen to be quite large. For 
the 7T°the theoretical value of the width is increased from r 7r o_ >2 7 = 7.74 cV to 

r>^ 27 = (8.25 ± .09) eV (2.28) 
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if we take fn+ffw = 1- Another correction to the chiral limit could originate 
from 7r°,?7 and 7/ mixing. It has been estimated in chiral perturbation theory 
[11] and would increase the decay rate by another 1.045 factor to 

r>^ 27 = (8.60 ± .10) eV (2.29) 

The current world average experimental value is 

r>^ 27 = (7.75 ± .60) eV (2.30) 

The large dispersion of the experimental results however suggests that the 
quoted error is underestimated [11]. The outcome of the PRIMEX experiment 
at Jefferson Lab [13], which aims at a precision of 1.5% in the measurement of 
the 7T° decay rate, is thus eagerly expected in order to clarify the situation. 

In order to extract information from our results on the 77 and rj more input 
has to be added 

An analysis of the 77, 77'mesons mass matrix together with chiral perturbation 
theory [10] yields for the mixing angle 9 = —19.5°, this value and the experi- 
mentally measured decay rates, when inserted in equations (2.22), (2.23) and 
(2.27) imply 

^8 = 1.29 ±.09 , F =4.70±.45 (2.31) 

Our result for F 8 is consistent with that obtained from chiral perturbation 
theory,F§ = 1.25 [10], the value obtained fori*o is quite larger than the ones 
appearing in the litterature.lt is worth noting that a pure gluon component 
present in 77, 77' would yield the value appearing in eq.(2.31) as an upper limit for 
F [A]. 

It is unfortunate that F$ and Fq are not directly linked to any physically 
measurable quantity. Data exist however on the transitions 77* — ► P [12]. 
We shall examine in the next section the implications of our results on the 
corresponding form factors . 

3 Form factors of the transitions 77* — ► rr , 77 ,7/ 

In the deep euclidean region we have [6] 

F p{ s,Q\t) = -2.{s-ml).j Q d XX wj o dy {y{x _ s + xQ2 f_ q2) _ xq2+m) 

(3.1) 

t = qf, Q 2 = —q\ and s = p 2 as before. In order to be able to compare with ex- 
periment, we need to evaluate F(m|,, Q 2 , 0). The method used in the preceding 
section,, modified to take into account the fact that we have now only a single 
pole in the vector meson channel now yields 

F(m 2 P ,Q 2 ,0) 
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(3.2) 

an expression valid for large Q 2 . Inserting expression (3.1) in the equation 
above , it takes some straightforward algebra in the complex plane to evalu- 
ate F(m 2 , Q 2 , 0). The deep inelastic limit Q 2 — > oo is readily obtained 

g.^(<g.,0)^<t^M±4pB (3.3) 

(m( - m 2 ).(m 2 - m 2 ) 

The corresponding expressions for the r\ and rf have to be multiplied by the 
factors 2 ^|™° ) and (2^» - 2v ^ OBfl ) respectively. Numerically 

then 

Q 2 .^o(m 2 ,Q 2 ,0) 

8^ (3 ' 4) 

g^^ ^ ( , 72 , ± , 4 ,^ + ^, ( , 6) 

From the experimental widths and the values of 6*, F , Fg given by equation 
(2.31) we get 

Q 2 .FJm 2 , Q 2 ,0) 

g^f j ^-72±-07 (3.7) 

Q 2 .F,,(m 2 Q 2 ,0) 



.98 ±.14 (3.8) 



The outputs of equations (3.4), (3.7) and (3.8) together with the experimental 
data [12] are shown on the figures 

4 Conclusion 

The chiral symmetry breaking corretions to the decays 7r°, r/, rf — > 27 have been 
estimated in a model independant way and found to be quite large. In the case 
of 7T° — ► 27 the theoretical value is onsidertably increased and this lends a great 
significance to the upcoming PRIMEX experiment [13]. 

In addition the form factors of the transitions 7*7 — + 7r°, 77, 77' were evaluated 
and shown to compare favourably with the available experimental data. 
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Figure Captions: 

Figl: The form factors of the transitions: 

a) 77* — > 7T° , b) 77* — > 77 , c) 77* — > 77' as obtained from eq. (3.2). 
The vertical lines represent the data as taken from [12] 
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